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ABSTRACT 

We present a pilot study on the origin and assembly history of the ICL for four galaxy clus¬ 
ters at 0.44 < z < 0.57 observed with the Hubble Space Telescope from the Cluster Lensing 
and Supernova Survey with Hubble (CLASH) sample. Using this sample of CLASH clusters 
we set an empirical limit on the amount of scatter in ICL surface brightness profiles of such 
clusters at z=0.5, a mean of 0.24 mag arcsec -2 for 10<r<110 kpc, and constrain the progeni¬ 
tor population and formation mechanism of the ICL by measuring the ICL surface brightness 
profile, the ICL colour and colour gradient, and the total ICL luminosity within the same 
radial range. This scatter is physical - it exceeds the observational errors, straightforward ex¬ 
pectations from the range of cluster masses in our sample, and predictions based on published 
evolutionary models for the variance attributable to the redshift span of our sample. We asso¬ 
ciate the additional scatter with differences in ICL assembly process, formation epoch, and/or 
ICL content. Using stellar population synthesis models we transform the observed colours to 
metallicity. For three of the four clusters we find clear negative gradients that, on average, 
decrease from super solar in the central regions of the Brightest Cluster Galaxy (BCG) to 
sub-solar in the ICL, under the assumption that the age of the intracluster stars is >11 Gyrs. 
Such negative colour (and equivalently, metallicity) gradients can arise from tidal stripping 
of L* galaxies and/or the disruption of dwarf galaxies, but not major mergers with the BCG. 
We also find that the ICL at 110 kpc has a colour comparable to m*+2 red sequence galaxies, 
suggesting that out to this radius the ICL is dominated by stars liberated from galaxies with 
L>0.2 L*. Finally, we find ICL luminosities of 4-8 L* between 10<r<110 kpc for these clus¬ 
ters. Neither dwarf disruption or major mergers with the BCG alone can explain this level of 
luminosity and remain consistent with either the observed evolution in the faint end slope of 
the luminosity function or predictions for the number of BCG major mergers since z=l. Taken 
together, the results of this pilot study are suggestive of a formation history for these clusters 
in which the ICL is built-up by the stripping of >0.2L* galaxies, and disfavour significant 
contribution to the ICL by dwarf disruption or major mergers with the BCG. 

Key words: galaxies: clusters: general, galaxies: elliptical and lenticular, cD, galaxies: evo¬ 
lution, galaxies: formation 


1 INTRODUCTION 


Galaxy clusters lie at the crossroads of astrophysics and cosmol¬ 
ogy, probing the influence of dark matter and dark energy in the 
evolution of large-scale structures in the universe. There exists a 
large body of research on cluster galaxies, the intracluster medium 
(ICM), and dark matter ( Kra vtsov & Bo rgani (2012), references 
therein). The properties of the intracluster starlight (ICL), the light 
from stars bound to the cluster potential but not to individual galax¬ 
ies, remain less well-determined than those of the other cluster 
components. In recent years, several deep ICL surveys of nearby 
clusters have bee n completed ([Gonzalez et al. |2000 [ Feldmeier 


|et al.||2002l|2004[ [Gonzalez et al.||2005||Krick et al.| |i006 Krick 


|& Bernstein) |2007l >. Individual stars in the ICL, such as intraclus¬ 
ter planetary nebulae (IPNe), supernovae (ISNe), and red giants 
(IRGs), have been used to constrain the kinematics of the ICL as 


well as the total luminosity of intracluster stars ( 

Feldmeier et al. 

1998 

Durrell et al. 

2002 1 Arnaboldi et al. 

2004, 

[Feldmeier et al. 

2004 

Gerhard et a 

l.| 2005; McGee & Bal 

ogh 2010, Sand et al. 

2011 

). Deep observations of Coma, Virgo, and other intermediate 


redshift clusters show distinct tidal tails and bridges, suggesting 
that the ICL grows during dynamical exchanges and can be used as 


a measure of the evolutionary stage of the cluster ( Feldmeier et al. 

2004 

Gerhard et al. 2005, Gonzalez et al. 2005, Mihos et al. 

2005 

Krick et al. 2006, Krick & Bernstein 2007, Rudick et al. 

2011 
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2 T. De Maio et al. 2014 


Because ICL formation is closely linked to the processes of 
cluster assembly, its study provides a means by which to under¬ 
stand the processes involved in the evolution of galaxies in clusters. 
Four scenarios have been put forward as possible mechanisms by 
which stars are added to the the ICL: (1) the disruption of dwarf 
galaxies as they fall into the cluster potential, (2) the tidal stripping 
of the outskirts of L* galaxies, (3) violent relaxation after major 
mergers between galaxies, including the central BCG, and (4) in 
situ star formation. Existing observational constraints on the envi¬ 
sioned formation scenarios are limited and theoretical models have 
yet to reach a consensus on which mechanism dominates the ICL 
formation process. 

The analytic models of [Purcell et al. |2007| and simulations 
of |Murante et al.| ( |2007| > and |Conroy et al.f i 2007]) posit that the 
majority of the ICL is built up either through mergers and/or inter¬ 
action with the BCG or through the shredding of dwarf satellites. 
Tidal stripping contributes only a small percentage of the ICL in the 
hydrodynamical simulations of Murante et al. ( 2007 ) and Sommer- 
Larsen et a l. ( 2005) whereas it is a more dominant mechanism in 
Rudick et~ak](|2009 1 and Watson & Conroy| ( [2013| )’s n-body simu¬ 
lations and |Laporte et al. ( |2013| )’s dark matter only simulations. In 
situ formation is less favoured as a means to produce intracluster 
stars. [Melnick et ak] ( |2012) ’s observations of the ICL of an inter¬ 
mediate redshift cluster attribute only 1% of the ICL mass to this 
population of stars despite the predictions ol jPuchwein et al.|(2010| 
that suggest ~30% of the intracluster stars are formed in situ. 

The study of the ICL is potentially significant in the con¬ 
text of resolving classic debates about the drivers of galaxy evo¬ 
lution in dense environments. Potentially important processes in¬ 
clude galaxy-galaxy tidal interactions and mergers, particularly in 
in-falling groups and subclusters (Zabludoff et al. (1996); Za blud-| 
off & Mulchaey ( 1998}, now called "pre-processing"), encounters 
with the cluster potential, galaxy harassment, and ram pressure 
stripping (e.g., |Park & Hwang] ( |2009| >; |Smith et~aT] ( |2010| ); |Wez-| 
gowi ec et al.| ( |2012| >). If the ICL is dominated by stars stripped 
from galaxies, it offers something unique: a signature of the his¬ 
tory of galaxy tidal interactions. By measuring the properties of the 
ICL over a range of redshifts and cluster masses, it is possible to 
constrain the importance of galaxy interactions in cluster evolution. 

In this work we use the surface brightness, infrared colour, and 
metallicity profiles of four massive galaxy clusters at intermediate 
redshift to comment on the progenitor population of the ICL and 
the dominant formation mechanism of the ICL. In Qwe describe 
the cluster sample used for this study and in ^3] we present the data 
reduction process, including treatment of the sky background. In ^4] 
we derive the ICL surface brightness and colour profiles. We next 
present the metallicity gradients in ^5]and discuss characteristics of 
individual clusters in ^6]before presenting the physical implications 
and conclusions of our results in ^7] and ^8] Throughout the paper 
we use a cosmology with H 0 =71 km s -1 Mpc -1 , £Im= 0.27, A=0.73. 


2 THE SAMPLE 

Our sample consists of clusters from the Cluster Lensing And 
Supernova survey with Hubble (CLASH, |Postman et al.| ( [20T2] )). 
An HST Multi-Cycle Treasury program, the CLASH sample con¬ 
sists of 25 galaxy clusters imaged in 16 different filters with 
WFC3/UVIS, WFC3/IR and ACS/WFC from 2000 A to 17,000 A. 
The 524 orbits of observations for this program began in Cycle 18, 
with the program scheduled for completion during Cycle 20. 

The CLASH sample is drawn mainly from the Abell and 


MACS cluster cat alogs ( [AbeTH [1958) |Abell et al.| [1989) |Ebeling| 
et al. |2007 2010| ) and contains clusters with masses of 5xl0 14 - 
3xl0 15 M 0 that range in redshift from z = 0.18 to z = 0.90 with 
a median redshift of z = 0.4. We refer the reader to postman et al.| 
P012| ) for additional details. In this paper, we focus on a subsample 
of four of the first observed CLASH clusters. With this subsample 
(Table [TJ we will study the evolution of the ICL for the upper mass 
range of galaxy clusters. 

While ICL analysis is not one of the main science goals of 
CLASH, the data set is particularly well-suited for the task. Intr¬ 
acluster light is intrinsically a low surface brightness component 
of galaxy clusters and thus benefits greatly from HST’s low sky 
background relative to ground-based observatories and instrument 
stability. We focus on the IR observations from CLASH because 
the stellar populations of the ICL are old and red and the spectral 
energy distributions are significantly redshifted. Utilizing infrared 
filters allows us to probe the ICL out to lower surface brightnesses 
and greater cluster radii than ever before, which enables us to view 
the ICL outside of the region dominated by the cluster BCG. 

Because sky subtraction is a critical step in the data reduction 
process, one aspect of CLASH that is not ideal for ICL analysis is 
the small field of view of WFC3. To minimize contamination by 
the ICL in the background estimation, we exclude the inner 300 
kpc of the cluster, which limits the amount of sky area available for 
background measurement. We develop a sky subtraction method 
that produces robust background estimations, despite the small field 
of view of HST, which we describe below in jj |3.1.3| 


3 REDUCTION AND ANALYSIS 

Using the standard calwf 3 package with reference darks, biases, 
and flats, we performed basic calibration for the data reduction. 
The first step of processing the data is performed as described in 
Schrabback et al. (2010). We use MultiDrizzle, which uses time- 
dependent distortion solution from Anderson (2007} to produce 
distortion corrected exposures. Each exposure is then interactively 
aligned by cross-correlating the positions of compact sources and 
applying residual shifts and rotations to the exposures and fed back 
into MultiDrizzle. The default values for MultiDrizzle’s cosmic ray 
rejection are adjusted slightly to allow for PSF variation due to tele¬ 
scope breathing, which can cause central stellar pixels to be flagged 
as cosmic rays. However, we turn off sky subtraction to avoid re¬ 
moving the ICL and instead rely upon our own method to apply this 
vital step (see j |3.1.3) . We use Scamp and SWaRP |Bertin[[2010a|b} 
for a higher precision WCS alignment between epochs and differ¬ 
ent filters. A standard tangential projection is used in Scamp and 
LANCZ0S3 resampling is employed by SWaRP. A planar fit en¬ 
sures minimal distortions from these two processes. Neglecting this 
step results in a ~2-4 pixel misalignment. For our highest-redshift 
cluster this error would correspond to a misalignment of 2.6 kpc 
and introduce artificial features in the inner colour profile of the 
ICL. 


3.1 Key Considerations 

Measurements of low surface brightness features like the ICL are 
intrinsically limited by background fluctuations. The two most 
dominant sources of uncertainties are flat fielding and background 
subtraction. Below we describe in detail our methodology for min¬ 
imizing uncertainties and biases associated with each. 
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Table 1 . Cluster Properties. 


Cluster 

Z c/«5 

M 5 ooc (10 14 M 0 ) 

M 2 ooc (10 14 M g ) 

Source 

MACS 1206.2-0847 

0.440 

10.6 + 2.1 

15.9 + 3.6 

Umetsu et al. 

|2014| 

MACS0329.7-0211 

0.450 

7.7 + 1.1 

10.0+1.5 

Umetsu et al. 

2014 

MACS 1149.6+2232 

0.544 

14.2 + 3.4 

25.4 + 5.2 

Umetsu et al. 

2014 

MACS2129.4-0741 

0.570 

10.6+1.4 

- 

Mantz et al. i 

20101 


3.1.1 Flat-Fielding 

Observations and analysis should be designed such that the error in 
photometry due to flat-fielding is negligible compared to the uncer¬ 
tainty due to sky subtraction. Super sky flats, high precision flat 
fields produced with stacks of images first masked of all astro¬ 
nomical objects, are used as the pipeline flats in the reduction of 
our sample ( |Pirzkal et al.[ |2011| . Across the whole WFC3 detec¬ 
tor these flats produce an RMS error in the photometry of 0.7% 
and a peak-to-peak uncertainty of 2.0+1.9%. In the central 700 
pix 2 of the CCD (excluding a border 1/8 the size of the detector) 
the RMS is 0.5%, with a peak-to-peak uncertainty of 1.5+1.6% 
( |Dressel[ |2011| ). For high surface brightness studies this level of 
systematic variation is insignificant; however, for the ICL this level 
of variation can introduce significant error. 

Two facts reduce the effects of low-frequency features in the 
flat field in our analysis. First, the CLASH field observations were 
taken at a variety of position angles. Field rotation distributes the 
observations over more pixels, thus averaging over large scale vari¬ 
ation. Second, we focus on radially averaged quantities centered on 
the BCG (see g4.2| ) so within a bin the error in photometry due to 
the field-fielding drops by a factor of 1 / ^/N, where N is the number 
of uncorrelated pixels in each azimuthal bin. 

To verify that the uncertainty in photometry due to flat-fielding 
is not a dominant source of systematic uncertainty, we create pre¬ 
liminary ‘delta’ flats for each filter by stacking WFC3/IR obser¬ 
vations bracketing the observation dates of these CLASH clusters. 
We drizzle these delta flats in a manner identical to the drizzling of 
the science images (same number of input images, with the same 
rotation). From these drizzled flats we compute for each filter a 
radially averaged profile of the residual flatness variations using 
the same methodology as described in ^4] We take the maximum 
amplitude of the radial flatness variation in the delta flats as the 
level of uncertainty introduced to the photometry due to variations 
in the flatness. Multiplying this variation by the measured back¬ 
ground level for each epoch of data, we find the variation of the 
measured surface brightness due to flat fielding uncertainty. For all 
clusters in both FI 10W and F160W, the systematic uncertainty due 
to flat-fielding is sub-dominant to the uncertainty in the measured 
background level. 

3.1.2 Masking 

To accurately measure the ICL we must avoid contamination. The 
extended light from galaxies and stars smoothly merges into the 
light from the sky, so separating the ICL from these sources is a 
non-trivial challenge. The most robust approach is to mask out all 
other sources, making sure not to too aggressively mask sources 
and leave too few pixels for ICL measurements, or inadequately 
mask sources and leave contamination from stars or galaxies in 
with the ICL. 

We approach this challenge via the following process. First, 
the sources are initially masked using ellipses extended out to 5 
times the semi-major and semi-minor axes determined by Source 


Extractor (Bertin & Arnouts 1996). A minimum detection area of 
five pixels and 2cr detection limit per pixel are required. Near the 
BCG we do not use the automated Source Extractor masking. In¬ 
stead, we manually mask any sources within 10" of the BCG by 
estimating the extent of the object by-eye and conservatively ex¬ 
tending the mask radius to ensure no galaxy light remains in the 
final ICL sample. To ensure that the extended wings of stellar pro¬ 
files are sufficiently masked, we enlarge stellar masks to 4" or to the 
extent of their diffraction spikes (~ 10-12" for the brightest stars). 
For even the brightest stars, this approach ensures that the light in 
the PSF wings is less than 2% of the sky level at 10", which is be¬ 
low the uncertainty in the measured background level and thus does 
not contribute significantly to the measured ICL signal. Each filter 
image is masked and the final cluster mask is produced by combin¬ 
ing the individual filter masks, thus the same masks are used for 
both filters of a given cluster, enabling self-consistent colour mea¬ 
surements. This method of masking is similar to those of |Jee 1 2010) 
and |Krick et al.| ( [2006) . See Figures [T|4| for pre- and post-masking 
examples. 


3.1.3 Background Subtraction 

While HST lacks the atmospheric emission that hinders ground- 
based observations, there is still background emission that must 
be considered. For pointings near the bright Earth limb, earthshine 
dominates over the contamination from zodiacal light and can pro¬ 
duce a significant large-scale gradient in the background levels. It 
is also possible for HST’s optical system to scatter stray light from 
bright sources outside the field of view into images ( [Dresse H |2ott) . 
These effects cause low-level variation in the background among 
epochs. Because of this variation we cannot use the available paral¬ 
lel observations in F125W and F160W to estimate the background 
in our science images as they were not taken synchronously with 
the science observations. Instead, we determine the background for 
each epoch of data individually from the target image. 

We first fit a plane to the sky, subtract off the gradients in x 
and y, and leave the normalization to be fit more rigorously for 
each epoch. This plane is fit using the unmasked background re¬ 
gion beyond 300 kpc from the BCG. Next the image is sliced into 
15 degree wedges, again excluding the inner 300 kpc and masked 
regions. For each wedge we determine the 3cr clipped median of 
the unmasked pixels in that wedge using the standard deviation of 
all sky values within a wedge as the lcr error values. The weighted 
average of all slices of a given epoch is then the final value used for 
the background level of that epoch. 

We find that the colour profiles produced by first subtracting 
the fitted background gradient are consistent with those produced 
without first subtracting the fitted planar background to lcr, sug¬ 
gesting that spatial gradients are not a significant driver of any mea¬ 
sured large-scale gradients. We also look at the time-dependent na¬ 
ture of the measured background gradients and find them to be sta¬ 
ble across the observations of a given cluster in a certain filter. We 
do this by differencing epochs of observations in a filter and fit- 
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ting the residual background gradient using the same least-square 
technique. The residual gradient results in a total edge-to-edge dif¬ 
ference in signal that is below the level of the measured background 
uncertainty. 

The wedge method of determining the background has the ad¬ 
vantage that any residual background gradient is averaged out by 
sampling the sky symmetrically. Because the sky is time-variable, 
we keep the data grouped by observation date in our analysis rather 
than combining after background subtraction. This choice has the 
additional benefit of allowing us to constrain the systematic errors 
between filters and observation dates. See Table[2]for the measured 
sky values with uncertainties. 

To test the robustness of our results, we use GALFIT ( |Peng| 
|et al.[ |2002| to model the background. We first mask out the 
BCG+ICL component to 30” on a masked, smoothed image, allow 
GALFIT to model the background with a planar model, and then 
renormalize the background values to have a median value equal to 
that found from the manual method described above. This alternate 
treatment yields results that are consistent at <1 cr to results from 
the manual method of sky subtraction. To use as few steps as pos¬ 
sible for the reduction process, we use the manual method of sky 
subtraction only. 


4 SURFACE BRIGHTNESS AND COLOUR PROFILES 
4.1 Surface Brightness Profiles 

Because the flat-fielding precision precludes a 2-D analysis of the 
surface brightness distribution at large radii, we produce radially 
averaged surface brightness profiles in FI 10W and F160W for each 
cluster, maintaining a roughly uniform signal to noise by measuring 
the flux in logarithmically spaced annuli. We take the median value 
of radii for unmasked pixels in a given bin as the bin location. This 
can lead to slight differences (~Alog(r[kpc])=0.007 in the locations 
of the bins between different clusters. Any masked pixels are ig¬ 
nored in flux calculations and are not replaced. In the top panel of 
Figure [6] we present the raw observed F160W surface brightness 
profiles, terminating them when the error in the surface brightness 
in F160W is >0.2 mag arcsec -2 or r<300 kpc. 

Extending the background apertures inward to 300 kpc while 
simultaneously measuring the surface brightness profiles out to 300 
kpc has the potential to result in over subtraction of the profile. 
However, we find that the impact of including the region at 300 kpc 
in the background annulus is negligible. This is mainly because 
most of the weight in determining the sky level comes from larger 
radii, as the number of pixels contributing to the background es¬ 
timate goes as r 2 . We have compared those background values as 
measured from 400 kpc outward to those from 300 kpc outward 
and find that they are the same within lcr. 

In the central panel of Figure [6] we correct for cosmological 
dimming, enabling a more direct comparison of the profiles. Fi¬ 
nally, in the bottom panel we correct the observed surface bright¬ 
ness profile not only for cosmological dimming but also for evolu¬ 
tion and passband (e+k) shifts. We use EZGAL ( jMancone ^Gon¬ 
zalez 2012) to produce e+k corrections for a Bruzual & Chariot 
(2003) (hereafter BC03) model with solar metallicity at a formation 
redshift of three (z/=3) and Chabrier IMF; see Table[5]for exact pa¬ 
rameters selected for each model. Using the evolution and passband 
corrections from the other available models has negligible effect on 
the final surface brightness profiles. Figure[7]is the same as Figure 
[6]but for the FI 10W filter. 


In Figure [8] we present the RMS scatter of the surface bright¬ 
ness profiles in bins of dlog(r[kpc])=0.15 in F160W (top) and 
FI 10W (bottom) for each successive correction (observed profiles 
in blue, //-dimming corrected in red, and //-dimming, passband, and 
evolution corrected as circles in black). At r >10 kpc, where the 
ICL is expected to dominate the surface brightness profiles, the 
scatter remains between 0.20-0.37 mag arcsec -2 in both filters, with 
a mean value of 0.27 (0.29) mag arcsec -2 in F160W (F110W). A 
single representative error bar for the uncertainty in the RMS is 
included on the right-most point for both filters. To factor out the 
differing masses of each cluster we scale the surface brightness pro¬ 
files of all clusters to a uniform mass of 1 x10 15 M o and find that the 
total RMS scatter in the profiles does not decrease. This suggests 
that other factors, such as the specifics of an individual cluster’s 
assembly history and ICL content, affect the scatter of observed 
surface brightness profiles. 


4.2 Colour Profile Extraction 


We present the BCG+ICL colour profiles that we produce by sub¬ 
tracting the surface profiles of each of the two filters bin by bin in 
Figure [9] The same masks are used for each filter, thus any result¬ 
ing features in the colour profile are not due to a lack of symmetry 
in masking. We plot all data to F160W=26 mag arcsec -2 , which 
effectively limits our profiles to r<110 kpc. Beyond this point the 
colour profiles become too noisy for use in robustly determining 
ICL properties. 

Due to the sensitivity of ICL measurements to background 
subtraction and the time variability of the sky, our largest source 
of systematic uncertainty is from sky brightness variation among 
epochs. The background level can vary by as much as 53% in 
F160W and 40% in F110W between observation epochs of a clus¬ 
ter in a given filter. (See Tablej3]for details of sky level variation 
per cluster per filter and Table[2]for background surface brightness 
levels with uncertainties for each epoch of data.) 

Large-scale background structure is evident across all epochs 
of data for both MACS 1149.6+2223 and MACS 1206.2-0847, 
which is also in a region of high zodiacal light. Additionally, 
there is known Galactic cirrus north-east of the cluster centre of 
MACS2129.4-0741 ( [Postman et al.[ |2012| >. The increased uncer¬ 
tainty in the colour profiles of these clusters at radii > 50 - 100 
kpc is a reflection of the presence of this structure. While both 
MACS 1206.2-0847 and MACS 1149.6+2223 show large-scale gra¬ 
dients in the background level, these gradients are stable across all 
epochs of data and thus do not alter the observed colour gradient 
beyond any effect that remains after symmetrically sampling the 
sky to determine each epoch’s background level (as described in 
g3. 1.3 1 above). 

To understand the effect of systematic errors between different 
filters and observation dates we analyse each epoch of data individ¬ 
ually and produce multiple colour profiles per cluster. We then av¬ 
erage together all the profiles produced with a single FI 10W epoch 
of data, e.g. 


yA (FllQWi - F16QWj) 
<F\\m-F\6m> i =2_ l — 1 


N 

N 


F\ 1 0W, - ^ 


F 160 Wj 

N 


( 1 ) 


We use the standard deviation of the mean as an estimate of the un¬ 
certainty because our systematics far outweigh our statistical errors. 
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Figure 1 . Single WFC3 F160W image of MACS 1206.2-0847. The circle is centered on the cluster BCG and has a radius of 200 kpc. The ICL cannot be 
traced beyond this radius because it is too faint compared to the background signal. Black regions have surface brightness brighter than 25.7 mag arcsec -2 
and white contours are overlaid in the central region at 21-24 mag arcsec -2 levels in steps of 1 mag arcsec -2 . The right image is the same as at the left but 
after masking all galaxies and stars. Masked regions are excluded and the average of the sigma-clipped median of the remaining pixels in each background 
region is the background value for the image. The -400 kpc low surface brightness bridge from the BCG to the north-east is a distinct feature, shown with a 
27 mag arcsec -2 level contour in black. Large-scale background structure is evident across all epochs of data for MACS 1206.2-0847; in the bottom corners of 
the above images, excess luminosity in the background is visible. 



A tabular version of the final FI 10W-F160W profile with lcr error 
values in dlog(r[kpc])=0.15 binning of each cluster is presented in 
Tabled] 

We investigate the ICL colour gradient by creating a single, 
averaged profile from all the individual colour profiles of each clus¬ 


ter. To maintain a uniform radial range for all clusters and produce 
a robust analysis, we fit only those radii for all clusters where the 
F160W surface brightness is < 26 mag arcsec -2 , which restricts our 
fits to r<110 kpc. Further, we do not fit inside 10 kpc, where the 
BCG dominates the surface brightness profiles. Using a weighted 
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Figure 3. Same as Figure[l]but for MACS 1149.6+2223. Black regions have surface brightness brighter than 26.4 mag arcsec 



Figure 4. Same as Figure[l]but for MACS2129.4-0741. Black regions have surface brightness brighter than 25.7 mag arcsec 2 . 


linear fit to constrain each profile for 10 < r [kpc] < 110, we find 
that the measured colour gradients of three of the four clusters are 
negative at the >3cr level (Table[6]). To derive these values we use 
radial profiles with dlog(r[kpc])=0.01 in order to avoid ambiguity 
in the discretization of the bins which may drive the fitted val¬ 
ues. In Figure [lO] we plot the observed colour gradients in order 
of increasing slope, top to bottom, and overlay the best fit, offset in 
colour for clarity. The large points are the colours of each cluster 
in dlog(r[kpc])=0.15, and the fitted lines are the results from fit¬ 
ting the dlog(r[kpc])=0.01 profiles. All clusters, with the exception 


of MACS0329.7-0211, trend toward significantly bluer colour at 
larger cluster radii. To test the time stability of the colour gradients 
we also fit the colour profiles resulting from averaging together the 
F110W-F160W profiles generated with each epoch of F110W and 
compare the resulting slope to that produced from the average of 
the colour profiles from all of epochs. The resulting slopes from all 
fits are consistent for all clusters, indicating that the reported colour 
gradients are not time dependent. 

We next apply passband and evolution corrections, derived 
from a BC03 model for a solar metallicity stellar population formed 
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Figure 5. The inner 200 kpc of MACS 1206.2-0847, MACS0329.7-0211, MACS 1149.6+2223, and MACS2129.4-0741 (left to right) in one epoch of F160W. 
White contours are at 21-24 mag arcsec -2 in 1 mag arcsec -2 increments and the black circle is at 200 kpc. Black regions are brighter than 25.7, 26.7, 26.4, 
and 25.7 mag arcsec -2 for each cluster, left to right, respectively. 


Table 2. Division of Observations. 


Cluster 

z 

Filter 

Date 

PA 1 

[deg] 

Exposure 

[sec] 

Sky 

mag arcsec -2 

SSky 

mag arcsec -2 

MACS 1206.2-0847 

0.44 

F160W 

4/3/11 

-44.33 

1005 

23.11 

28.65 




6/21/11 

-44.33 

1508 

22.85 

28.85 




5/2/11 

-15.33 

1005 

23.27 

28.87 




7/5/11 

-20.33 

1508 

22.69 

28.62 



FI 10W 

4/3/11 

-44.33 

1508 

23.35 

29.06 




6/21/11 

-44.33 

1005 

23.09 

29.37 

MACS0329.7-0211 

0.45 

F160W 

8/18/11 

133.67 

1005 

23.09 

29.66 




10/16/11 

133.67 

1508 

23.37 

29.35 




9/20/11 

155.67 

1005 

23.32 

29.37 




11/1/11 

155.67 

1508 

23.23 

29.15 



FI 10W 

8/18/11 

135.67 

1508 

23.27 

30.13 




10/16/11 

135.67 

502 

23.50 

30.45 

MACS 1149.6+2223 

0.544 

F160W 

1/16/11 

124.67 

1005 

21.79 

27.58 




2/27/11 

124.67 

1508 

21.79 

27.71 




12/4/10 

156.67 

1005 

21.46 

27.62 




3/9/11 

96.67 

1508 

21.79 

27.67 



FI 10W 

1/16/11 

124.67 

1508 

21.98 

28.52 




2/27/11 

124.67 

402 

21.81 

27.96 

MACS2129.4-0741 

0.57 

F160W 

5/15/11 

94.67 

1005 

21.38 

28.45 




6/25/11 

123.67 

1205 

21.81 

28.37 




6/3/11 

94.67 

1408 

21.65 

28.38 




7/20/11 

94.67 

1408 

21.76 

28.50 



FI 10W 

5/15/11 

94.67 

1408 

21.50 

28.52 




7/20/11 

94.67 

502 

21.94 

28.70 


1 Position angle from the reference aperture center. 


Table 3. Sky variation between epochs. 


Cluster 

Filter 

Sky variation 

% 

dSky 

mag arcsec -2 

MACS 1206.2-0847 

F160W 

53.1 

27.98 


FI 10W 

23.9 

28.81 

MACS0329.7-0211 

F160W 

26.0 

28.60 


FI 10W 

21.6 

29.89 

MACS 1149.6+2223 

F160W 

32.3 

26.89 


FI 10W 

15.2 

27.73 

MACS2129.4-0741 

F160W 

40.6 

27.67 


FI 10W 

39.7 

28.23 


at Zy—3, to the observed F110W-F160W colour profiles to bring 
them to their z=0 colours. The resulting profiles are shown in Fig¬ 
ure m where we plot each cluster slightly offset in radius for clar¬ 


ity. The mean colour of all four clusters in each dlog(r[kpc])=0.15 
bin is plotted in large black triangles. We perform a linear fit to 
the mean colour profile, resulting in a colour gradient of d(Fl 10W- 
F160W)(d log(r)) -1 =-0.055+0.010 mag arcsec -2 log(kpc) -1 . 


5 METALLICITY GRADIENTS 

Evolutionary models ( [Melnick et al.| |2012[ IRudick et al.[ |2006[ 


Sommer-Larsen et al. 1 12005| |Murante etltL[ |2004[ |2007| |Conroy 


|2007| > and observations ( |Feldmeier et al.||2004[|Krick et al. 


et al. 


20061 suggest that the formation of the ICL is closely related to the 


dynamical history of the cluster, thus its study provides valuable in¬ 
sight into the cluster formation process. One way to distinguish be¬ 
tween the many ICL and cluster assembly scenarios is to look at the 
metallicity of the ICL. The existence of a mass-metallicity relation 
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Table 4. Colour Profiles 


MACS 1206.2-0847 MACS0329.7-0211 MACS 1149.6+2223 MACS2129.4-0741 


log(r[kpc]) 

F110W-F160W 

log(r[kpc] 

F110W-F160W 

log(r[kpc] 

F110W-F160W 

log(r[kpc] 

F110W-F160W 

0.08 

0.489 + 0.031 

0.10 

0.496 + 0.018 

0.11 

0.497 + 0.038 

0.08 

0.468 + 0.026 

0.24 

0.484 + 0.016 

0.24 

0.442 + 0.059 

0.21 

0.494 + 0.020 

0.23 

0.469 + 0.024 

0.39 

0.523 + 0.036 

0.37 

0.438 + 0.016 

0.39 

0.531 +0.027 

0.38 

0.513+0.028 

0.53 

0.474 + 0.015 

0.53 

0.427 + 0.018 

0.56 

0.484 + 0.023 

0.54 

0.534 + 0.017 

0.69 

0.479 + 0.016 

0.68 

0.424 + 0.040 

0.70 

0.507 + 0.018 

0.69 

0.525 + 0.009 

0.84 

0.485 + 0.007 

0.84 

0.434 + 0.005 

0.84 

0.501 + 0.006 

0.83 

0.510 + 0.012 

0.98 

0.492 + 0.008 

0.99 

0.425 + 0.007 

0.99 

0.474 + 0.007 

0.98 

0.522 + 0.013 

1.14 

0.492 + 0.010 

1.14 

0.438 + 0.008 

1.14 

0.470 + 0.002 

1.14 

0.495 +0.010 

1.29 

0.491 +0.007 

1.29 

0.439 + 0.003 

1.29 

0.458 + 0.010 

1.29 

0.494 + 0.012 

1.44 

0.467 + 0.012 

1.44 

0.449 + 0.012 

1.42 

0.450 + 0.012 

1.44 

0.445 +0.012 

1.58 

0.452 + 0.006 

1.59 

0.437 + 0.008 

1.60 

0.456 + 0.015 

1.58 

0.477 + 0.038 

1.74 

0.423 + 0.010 

1.73 

0.439 + 0.006 

1.73 

0.403 + 0.024 

1.74 

0.469 + 0.064 

1.89 

0.422 + 0.020 

1.89 

0.441 + 0.015 

1.89 

0.413 + 0.056 

1.89 

0.471 +0.066 

2.04 

0.457 + 0.059 

2.04 

0.399 + 0.011 

2.04 

0.351 +0.109 

2.04 

0.460 + 0.166 


( [Skillman et al.[[T996} in cluster galaxies implies that one can iden¬ 
tify the characteristic mass of the progenitor population using the 
metallicity of the ICL. Intracluster stellar metallicities have been 
spectroscopically measured in Virgo and Fornax from several types 
of intracluster populations including RGB stars, planetary nebu¬ 
lae, and supemovae ([Durrell et al.l |2002| |Feldmeier et al.[ |1998| 
[Williams et al.[|2007]|Coccato et al.[|2011) . In general, these studies 

find that intracluster stars are between 2-13 Gyr old and have sub¬ 
solar metallicities (-0.8<[Fe/H]<-0.2) ( [Durrell et al.[|2002||Coccato| 
[eFaLl[20TTl >. We use the colour gradients of the integrated stellar 
population presented in j }4.2| to constrain the metallicity profiles of 
the ICL, modulo uncertainty in the age of the stellar population. 


5.1 Stellar Population Models 

To transform the observed colour gradients into metallicity profiles, 
we employ stellar population synthesis (SPS) modeling. We use 
five widely-used model sets, taking the models of |Bruzual & Char-] 
[lot] ([2003) as the fiducial to examine the effect of model choice on 
our results. These include an updated version of the BC03 model 
(often referred to as CB07) and the model set from |Maraston| ( |20()5] ) 
(hereafter M05). Both include detailed treatment of thermally pul¬ 
sating asymptotic giant branch stars (TP-AGB), which can have a 
potentially significant effect on the infrared emission of interme¬ 
diate age stellar populations. Additionally, we use Flexible Stellar 
Population Synthesis models from Conroy et al. ( 2009]) andfConroy 
|& Gunn| ( |2010| ) (hereafter C09), which treat some SPS inputs, such 
as IMF and stellar evolution uncertainties, as free parameters. We 
also use the models from |Percival et al.] ( |2009] > (hereafter BaSTI), 
which have a larger range of metallicities and cr-enhanced models. 
Table [5] lists available metallicities and IMFs for each model set. 
We primarily use a Chabrier IMF jChabrier) |2003|> for the BC03 
and CB07 model sets and Kroupa IMF ( Kroupa] |2001| ) for the C09, 


M05, and BaSTI model sets. We employ a simple stellar population 
for our star formation history. From these SPS models we derive 
passband and evolution corrections (see Q, as well as metallici¬ 
ties of the ICL given an observed colour (see m 

The clear blue-ward radial trend in the FI 10W-F160W colour 
profiles seen in Figure [lO] can be interpreted as either an age-radius 
or a metallicity-radius relation, however there exist both theoret¬ 
ical and observational reasons to expect that the dominant factor 
is metallicity. Observationally, local studies generally indicate that 
the intracluster stellar population is an old stellar population (>10 


the intracluster stars form from z~2 ( [Murante et al.[ |2004[ ) to z~3 
(Sommer-Larsen et al. 2005, Puchwein et al. 2010). Because of 
the general consensus that the ICL is an old stellar population we 
explore the interpretation that our observed colour gradients are due 
to changes in metallicity, noting the potential impact of the age- 
metallicity degeneracy. 


5.2 Metallicity Profiles 

We transform our observed colour gradients to metallicity profiles 
with each of the SPS models, using EZGAL ( |Mancone & Gonza^] 
|lez[ |2012| ), which outputs observable parameters from input SPS 
models. For a given model set, cluster redshift, assumed formation 
redshift, and star formation history, the metallicity of a given stel¬ 
lar population has a unique corresponding colour. Over the range 
of metallicities available for each model we build up a metallicity - 
colour function for each model set. The ICL metallicity gradient of 
a cluster is then found by interpolating the observed colour in each 
radial bin to a corresponding model-derived metallicity. 

The colour gradients in Figure [TT] directly translate into corre¬ 
sponding metallicity gradients under the assumption of a fixed age 
for the ICL. Using a formation redshift of z f =3 as a baseline, the 
resulting metallicity gradients are shown in Figure [12] We use the 
F110W-F160W colour to interpolate to metallicity as F110W is 
less noisy than F125W and FI 10W-F160W provides greater wave¬ 
length leverage. Of the five models, we show only M05, BC03, and 
CB07 in Figure [i~2| BC03 is our fiducial model and generally lies 
between the more extreme predictions of M05 and CB07. For clar¬ 
ity, we do not show the metallicity gradients for C09 and BaSTI, 
which fall between the models shown and exhibit the same general 
shape. The ar-enhanced BaSTI models produce gradients consistent 
with the zero ar-enhancement BaSTI models. We investigate the ef¬ 
fects of IMF choice on the produced metallicity gradients and find 
little effect; the gradients produced with a Salpeter IMF are consis¬ 
tent with those produced using a Chabrier or Kroupa IMF. 

We fit the metallicity profile of each cluster with a weighted 
linear fit between 10<r<110 kpc and find that the metallicity pro¬ 
file of MACS0329.7-0211 is consistent with being flat, while the 
other three clusters have a range of metallicity gradients of - 
0.510<d(Z/Z o )(d log(r)) -1 <-0.659. In Figure 13 we overlay the 


weighted linear fits on the interpolated metallicity points for a 
BC03 model with zy=3. Using the results from the BC03 model 
as a fiducial, we see that at radii <10 kpc all clusters have abun- 


Gyrs, 

Williams et al. 

(2007,; 

|Coccato et al. ( 2011); Melnick et al. dances >Z Q and all but MACS0329.7-0211 dip to solar abundance 

12012 

); Loubser et a 

(2009 

)). Theoretically, models predict that and below at larger cluster radii. 
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Table 5. Model Set Parameters. Adapted from |Mancone & Gonzale/|f201 2| , Table 1. The primary IMF model used for each stellar population synthesis model 
in our analysis is bolded. 


Model 

Metallicity (Z/Z Q ) 

u-enhanced 

# Metallicities 

Salpeter 

Chabrier 

Kroupa 

BC03 

0.005-2.5 

No 

6 

Yes 

Yes 

No 

CB07 

0.005-2.5 

No 

6 

Yes 

Yes 

No 

BaSTI 

0.005-2 

Yes 

10 

No 

No 

Yes 

C09 

0.01-1.5 

No 

22 

Yes 

Yes 

Yes 

M05 

0.05-3.5 

No 

5 

Yes 

No 

Yes 



10 100 


r [kpc] 

Figure 6. F160W surface brightness profiles of BCG+ICL as a function of 
radius in dlog(r[kpc])=0.15 bins, out to the point where the error in F160W 
surface brightness for each cluster is <0.2 mag arcsec _2 and r<300 kpc. The 
top panel shows the observed profiles with no passband or evolution cor¬ 
rection. The middle figure takes the observed profiles and corrects for cos¬ 
mological dimming. Finally, the bottom panel accounts for evolution and 
passband (e+k) corrections, assuming a BC03 model, as well as cosmolog¬ 
ical dimming. 


At low redshifts [Loubser et al.| j2009) find spectroscopically 
measured BCG metallicities of ~2 Z©. Looking to our results, and 
using BC03 as a fiducial, we find that all four clusters have simi¬ 
lar BCG metallicity of -1-1.4 Z Q for the BC03 model. Referring 
to Figure [12] we see that the M05 model set predicts metallicities 
closest to the observed range of Loubser et al. ( 2009]) and[Udman 
|etaLl ( [20T2l ). 

While there are many free parameters for each model set, the 
most significant for our analysis is the formation redshift. Our base¬ 
line assumption of Zf=3 is in agreement with the model predictions 



10 100 


r [kpc] 

Figure 7. The same as in Figure[6] but for FI 10W. 


and observations |Sommer-Larsen et al. 

2005, Puchwein et al. 

2010 

Williams et al. 2007, Loubser et a] 

. 2009, Coccato et al. 

2011 

Melnick et al. 2012). We present the metallicity profile of 


MACS 1149.6+2223 as a function of formation redshift for a BC03 
model (Figurep4| to illustrate how the choice of formation redshift 
affects the final metallicity gradient. From z f =2 to z f =6 we see a 
~0.2 Z/Z© decrease in the inferred metallicity. However, the down¬ 
ward trend in metallicity with radius remains, regardless of which 
formation redshift is used. This suggests that, despite the uncer¬ 
tainty in the ICL age, the trends exhibited by our metallicity pro¬ 
files are robust and that only when considering absolute metallicity 
does model choice play a noticeable role in the ICL metallicity. 
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10 100 


r[kpc] 

Figure 8. The RMS scatter between the surface brightness profiles of all 
four clusters in F16W (top) and F110W (bottom) for each successive sur¬ 
face brightness correction. The scatter between observed profiles are traced 
by the blue line, the scatter between //-dimming corrected profiles by the 
red line, and completely corrected (//-dimming & k+e corrected) in circles 
with a connecting black line. For 10<r<110 kpc the scatter remains be¬ 
tween 0.20-0.37 mag arcsec -2 in both filters, with a mean of 0.27 (0.29) 
mag arcsec -2 in F160W (FI 10W). The error bar on the largest radius point 
reflects the estimated uncertainty in the plotted RMS values. 



r [kpc] 


10 100 
r[kpc] 


Figure 9. F110W-F160W observed colour profiles, with no passband or 
evolution correction applied, for all four clusters in order of increasing red- 
shift, top to bottom. We produce 2 averaged colour profiles, one for each 
blue filter epoch of data (solid and dashed line, respectively). All radii with 
^F160W < 26 mag arcsec -2 are plotted with error bars. ICL measurements 
are sensitive to background subtraction and thus our largest source of sys¬ 
tematic uncertainty is from sky brightness variation among epochs of data. 
Differences between filters and observation dates drive the larger error bars 
and differences in colour seen in the colour profiles at large radii. 


6 NOTES ON INDIVIDUAL CLUSTERS 


Two of our clusters, MACS0329.7-0211 and MACS 1206.2-0847, 
are included in the list of eight CLASH clusters with possible sub¬ 
structure (Post man et al.|[20l2] >. The colour and metallicity gra¬ 
dients of the ICL of MACS0329.7-0211 may be flattened due to 
light that is contributed from a second giant elliptical galaxy located 
40" north-west from the central BCG. To test this supposition we 
mask the entire half of the cluster containing the second giant ellip¬ 
tical and compare the colour profiles before and after the additional 
masking. We find the colour profiles to be consistent, suggesting 
that if the two galaxies have interacted and deposited stellar pop¬ 
ulations into the ICL, those additions have been well-mixed with 
the rest of the ICL. postman et al.| ( |2012| ) note that MACS0329.7- 
0211 was initially classified as dynamically relaxed by | Allen et al. 
(20081 and |Schmidt & Alien] ( |2007| > but that later |Maughan et al" 
2008) found evidence of substructure. The possible dynamical his¬ 
tory that this substructure suggests could have flattened the stellar 
population gradient of the BCG+ICL component and produced the 
observed nearly constant colour and metallicity of MACS0329.7- 
0211’sICL. 


The dynamical state of MACS 1206.2-0847 remains unclear. 
Many signs point to the cluster being dynamically relaxed: the op¬ 
tical position of the BCG and X-Ray peak are coincident ([Post¬ 


man et al.] 2012| >, it has a symmetric X-Ray morphology (Gilmour 


etal. 


|2009 l, the dark matter, BCG, and X-Ray centroids are aligned 


pJmetsu et al.[ |2012| |Biviano et al.[ |2013| > and show no signs of 
recent merging activity. However, there are multiple signatures in¬ 
dicative of recent merger activity. The se signposts inc lude a high 
velocity dispersion (>1500 km s -1 ) ( [Gilmour et al. 2009 


Bi- 


viano et aL]|2013|), high X-Ray luminosity and temperature (Ebel- 
ing et al.[|2010[ postman et al.]|2012| >. Additionally, MACS 1206.2- 
0847 displays a low surface brightness structure (down to ~27 mag 
arcsec -2 ) extending from the BCG roughly 400 kpc to the north¬ 
west (See Figure [T] presotto et al.| ( |2014] ); |Umetsu et al.| ( [2012| >). 
Such a distinct feature, which is aligned with the BCG position 
angle and is traced by in-falling galaxies presotto et al.[ |2014j ), 
suggests recent dynamic activity. 


7 DISCUSSION 

7.1 ICL Progenitor Populations 

Here we consider the possible mechanisms of formation of the ICL, 
describing how the colour and metallicity gradients can be used as 
diagnostics. 


7.1.1 Formation Mechanisms 

The exact mechanisms of ICL formation remain poorly con¬ 
strained. Four distinct mechanisms have been proposed to explain 
the formation of the ICL: 

(1) Shredding of dwarf galaxies. As low-density galaxies fall into 
the cluster, they can be tidally disrupted by the cluster poten¬ 
tial or by encounters with cluster member galaxies. Their stars 
become unbound, responding to the cluster potential, and form 
the intracluster starlight. The most metal-poor dwarfs, which 
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Figure 10. The observed F110W-F160W averaged colour profile of each 
cluster, arranged by increasing slope, top to bottom. The large points are 
the colours of each cluster in dlog(r[kpc])=0.15 with error bars representing 
the standard deviation of the mean colour. Clusters are offset in colour by 
the amount given at the right side of the figure and the best-fit linear equa¬ 
tions resulting from fitting the dlog(r[kpc])=0.01 profiles for each cluster are 
overlaid. All clusters except MACS0329.7-0211 trend toward significantly 
bluer colours at higher radius with a colour gradient d(F110W-F160W)(d 
log(r)) -1 between -0.063 and -0.101 mag log(kpc) -1 (see Table[6j. 


Figure 11. The F160W-F110W colour profiles after passband and evolu¬ 
tion correction with a BC03 model for a solar metallicity stellar population 
formed at zy=3. Each cluster is slightly offset in radius for clarity. The large 
black points are the mean in each dlog(r[kpc])=0.15 bin with the standard 
error in the mean as the lcr error bars. The black line represents the linear fit 
found for the mean colour vs. radius relationship. The RMS scatter between 
the colour profiles of the four clusters is shown in the lower panel, which is 
low and roughly uniform from 10<r<l 10 kpc. 


are also the lowest mass due to the mass-metallicity relation, 
are disrupted at larger cluster radii than more massive satellites, 
thus forming a gradient in colour and metallicity with cluster 
radius (|Rudick et al.||2009||Melnick et al.||2012||Conroy et al.| 
[2007). 

(2) Tidal stripping. The tidal stripping of L* galaxies from inter¬ 
actions with other cluster members or the cluster potential can 
contribute to the build-up of the ICL. Because L* galaxies have 
radial metallicity gradients (Zaritsky et al. 1994, La Barbera 
|et al.[|20T2| ), the tidal forces of the cluster potential will liber¬ 
ate progressively more tightly bound (redder) stars as a galaxy 
approaches the cluster centre. The stars added to the ICL will 
therefore be somewhat more metal-rich at smaller cluster radii. 
See iMelmck et al.| <20T2] ); |Rudick et all ( [20091 |2006l [2010] ); 
|Puchwein et al.| (|2010|>; |Murante et al.|(|2007|>; |Conroy et al.| 
12007) . 

(3) Violent relaxation during mergers. The ICL may be built up 
via liberation of stars from galaxies undergoing violent merg¬ 
ers with the centrally dominant cluster galaxy. It has been sug¬ 
gested that 30-50% of stars of the merging body become un- 


bound and end up in the ICL 

Murante et al. 20071 Conroy| 

et al.| 20071 Lidman et al. 2012 

). Further, major mergers serve 


to flatten any existing stellar population gradient ( White 1980 

Kobayashi 2004, Di Matteo et al. 

2009, La Barbera et al. 

2012, Eigenthaler & Zeilinger 2013 

). Thus, if major mergers 


are a dominant process of ICL formation, then the BCG+ICL 
colour and metallicity gradients are expected to be shallow or 
near-constant. The metallicity of the ICL will depend on the 


masses of the galaxies that interact with the BCG. The more 
often massive galaxies merge with the BCG, the more metal- 
rich the intracluster stellar population will be. 

(4) In situ formation. Finally, a portion of the intracluster stel¬ 
lar population may arise from in situ star formation ( |Puch-| 
|wein et al.||2010) , wherein clumps of gas stripped from clus¬ 
ter galaxies remain cool and compact enough to form new 
stars. While the impact of ongoing star formation on the colour 
and metallicity gradient of the ICL+BCG is unclear, the pro¬ 
cess would change the mean age of the ICL stellar popula¬ 
tion, which must be taken into account to successfully inter¬ 
pret any observed colour trends. Little observational evidence 
exists suggesting that in situ stars are a significant population 
in the ICL. For example, using synthetic models for the ICL, 
|Melnick et al.| ( |2012| find that only 1% of the ICL consists 
of in situ stars in their study of a z=0.29 cluster. The lack of 
Type la intracluster stellar supernovae in z~0.1 clusters ( |Sand| 
|et al.[ |2011| ) further suggests that in situ star formation is not 
the dominant mechanism of ICL formation. 


Current simulations present a range of predictions for the rela¬ 
tive importance of these mechanisms. Some suggest that disruption 
of satellites and tidal stripping (mechanisms (1) & (2) above) may 
contribute significantly to the ICL stellar population ([Purcell et al. 


2007 Stanghelfini et al.[| 2006||Rudick et al.[|2009||Contini et al. 


2013 Laporte et al.[|2013| >, whereas others predict that the major¬ 

ity of the ICL is built up from violent mergers of cluster mem¬ 


bers with the centrally dominant galaxy (mechanism (3), Murante 
|et al.|([2007] >; |Conroy et al.| ( |2007| )). The results of |Contini et al. 
P013| ) favour the scenario in which stellar stripping and disrup- 
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Figure 12. SPS models are employed to transform the observed F110W-F160W dlog(r[kpc])=0.15 colour profile for each cluster to metallicity as a function 
of radius. Errors are those errors from the colour profiles in Figure[To| transformed to metallicity in the same fashion. In general, there is considerable spread 
between the metallicity profiles produced with the five stellar population models, though the basic shape of each profile remains the same. We show only 
BC03, CB07, and M05 for clarity. 


tion of satellite galaxies are the primary mechanisms by which the 
ICL forms. [Contini et al~|<|2Ql 3|) use updated semi-analytical mod¬ 
els based upon De Lucia & Blaizot ( 2007) to specifically address 
how the ICL forms and which progenitor population contributes 
most to the ICL stellar population. They suggest that mergers con¬ 
tribute relatively little to the ICL, in contrast with the results of 
|Murante et al.|p007) , who argue that tidally-stripped stars repre¬ 
sent less than 10% of the total ICL. 

Using their tidal stripping and satellite disruption models, 
|Contini et al.| ( |2013} find that 26% of the ICL starlight results from 


galaxies with stellar masses in the range of 10 10-75 -10 1L25 M© and 
that 68% of ICL starlight comes from galaxies with M* >10 10 5 
M 0 . Intracluster stars from low-mass galaxies (M* < 10 9 Mq) con¬ 
tribute negligibly to the total ICL population. This, they argue, is a 
consequence of low-mass satellites spending more time at the clus¬ 
ter outskirts, where tidal stripping and disruption are less efficient. 
This picture fits with observed cases of mass segregation in clusters 
e.g. |Andreon| ( |2002] ); |Biviano et al.| ( [2002] ); |Presotto et al.|(2012| ). 

Using the observed colour and metallicity gradients of the ICL 
presented here, we comment on the likely importance of the vari- 
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Figure 13. The interpolated metallicity profiles of the BCG+ICL produced 
using the BC03 model set for zy=3, over-plotted with their best fit linear 
relation to illustrate the similarity between clusters. The large points are 
from the metallicity profiles produced with the dlog(r[kpc])=0.15 colour 
profiles and the best fit relation shown is the result of a weighted linear fit 
to the dlog(r[kpc])=0.01 profiles. All four clusters show ~Z© metallicities 
within 10 kpc and either remain at Z© or dip to sub-solar metallicities by 110 
kpc. MACS0329.7-0211 is the only cluster to show a flat metallicity profile, 
which is likely due in part to the cluster’s unrelaxed dynamical state and 
contribution of ICL from the massive elliptical galaxy located at a projected 
40" north-west from the central BCG. 

ous mechanisms of ICL formation and on the dominant intracluster 
stellar progenitor population. 


7.2 ICL Assembly History 

Over the redshift range of our sample (z=0.44-0.57) the evolution¬ 
ary models of |Contini et al.| ( [2013| ) predict that the ICL grows by 
-40%, increasing from -35 to 50% of its present day mass for their 
"disruption" and "tidal" models (See their Figure 6). If we assume 
that 1) our clusters are different versions of the same evolving clus¬ 
ter with an age spread corresponding to the redshift interval of our 
sample and 2) the ICL mass growth corresponds directly to a rise 
in surface brightness, then both of their models predict a scatter 
of 0.15 mag arcsec -2 among the surface brightness profiles. As we 
have neglected all sources of scatter other than cluster redshift, the 
model estimates are strict lower limits. 

The observed profiles of the four clusters in our sample ex¬ 
hibit a mean scatter of 0.27 mag arcsec -2 between 10<r<110 kpc. 
At these radii, the physical scatter dominates over observational 
uncertainties. Furthermore, as noted in Section [Li] a simple mass 
rescaling does not reduce the measured scatter. Therefore, for the 
range of cluster mass probed here, we find no dependence of the 
normalization of the ICL profile on cluster mass. 

We next estimate the contribution of observational errors in 
the surface brightness profiles to the scatter. Quantifying the er¬ 
ror in the scatter empirically from the systematic variations among 



Figure 14. Effect of formation redshift on the metallicity interpolated from 
the F110W-F160W colour of the BCG+ICL of MACS 1149.6+2223 using 
the BC03 model set. We produce interpolated metallicity profiles using 
Zf=2-6, omitting errors bar on all but the zy=3 line for clarity. Changing 
the formation redshift from z/=2 to z/=6 affects the metallicity by a max¬ 
imum of only ~0.2 dex and, most importantly, preserves downward trend 
in metallicity with radius. This suggests that, despite the uncertainty in the 
age of the ICL, the general trends in metallicity that we observe are robust. 

our profiles is difficult due to the small sample size of our current 
sample. With this caveat in mind, we nevertheless use our data to 
estimate an uncertainty in the RMS values of -0.05 mag arcsec -2 , 
plotting an error bar for reference in the rightmost point in Figure 8. 
Many points are more than lcr above the model RMS of 0.15 mag 
arcsec -2 , and all lie above it. We conclude that the measured scat¬ 
ter does indeed exceed what is predicted from simple mass growth, 
suggesting a wider range of assembly processes, epochs, and/or 
ICL content. We leave a more detailed discussion for future papers 
that will contain our larger, complete sample of clusters. 

None of the existing theoretical studies publish the scatter 
among surface brightness (or colour) profiles as a function of red¬ 
shift or radius. We advocate for direct comparisons between models 
and theoretical predictions for the scatter in ICL build-up among 
clusters, particularly as the observational work expands to larger 
samples spanning a wider redshift range. Our future work with the 
full set of CLASH clusters ( [Postman et al.[ |2012| > and a sample of 
galaxy groups (HST Program 12575) is intended to provide this 
type of benchmark, as it will use systems across a range of red- 
shifts with masses that reflect the expected hierarchical build-up of 
clusters. 


7.2 .1 Colour and Metallicity Gradients as Diagnostics 

The colour and metallicity of the ICL reflects the dominant progen¬ 
itor population, encoding the formation history of the cluster in the 
ICL ( [Zibetti etaT]|2005[|Willman et al.[|2004||Murante et al.[|2004l 

|Sommer-L arsen et al. 2005). Purcell et al. ( 2007) predicts that clus¬ 
ters with a halo mass greater than 10 13 5 M© will have an ICL of so- 
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Table 6. ICL Colour, Metallicity Gradients and Luminosity 


Cluster 

z 

d(F 110W-F160W)(d log(r))-‘ 
[mag arcsec -2 log(kpc) -1 ] 

d(Z/Zo)(d log(r))- 1 
[log(kpc) -1 ] 

L(r<10 kpc) 
[L*] 

L(r<110 kpc) 
[L*] 

L(10 < r(kpc) < 110) 
[L*] 

MACS 1206.2-0847 

0.44 

-0.091 ± 0.006 

-0.659 ± 0.035 

1.11 

5.76 

4.65 

MACS0329.7-0211 

0.45 

0.000 ± 0.004 

-0.025 ± 0.018 

1.43 

8.04 

6.61 

MACS 1149.6+2223 

0.54 

-0.101 ±0.014 

-0.537 ± 0.062 

1.75 

10.26 

8.51 

MACS2129.4-0741 

0.57 

-0.063 ± 0.014 

-0.510 + 0.105 

1.47 

6.85 

5.38 


lar metallicity and the simulations of |Sommer-Larsen et al.| ( [2005| 
show ICL colours consistent with a progenitor population of sub- 
L*, E, and SO galaxies. 

The colour gradients of the ICL in our clusters signifi¬ 
cantly trend toward bluer colour at larger cluster radii, except for 
MACS0329.7-0211, which has a flatter profile (See Figure [TU]). [Zn] 
|betti et alT| ( |2005) similarly find a blue-ward gradient in their anal¬ 
ysis of stacked SDSS clusters and |Montes & Truj iTTo] < |20 1 4| > show 
that the Frontier Fields cluster Abell 2744 has a bluening colour 
profile in both g-r and i-J. These observations strongly favor the 
formation mechanisms that produce distinct metallicity gradients: 
tidal stripping and disruption of the L* and dwarf galaxy popula¬ 
tions. 

In addition to the F110W-F160W colour gradients, we also 
use the absolute metallicities and metallicity gradients of the ICL 
in these clusters to constrain the progenitor population and forma¬ 
tion mechanism of the ICL (See Figure[l3]). To this end we perform 
weighted linear fits to the metallicity profiles from the BC03 model 
for 10<r<110 kpc. MACS 1206.2-0847, MACS 1149.6+2223, and 
MACS2129.4-0741 all exhibit steep metallicity gradients of, on 
average, d(Z/Z Q )(d log(r)) -1 - -0.6 log(kpc) -1 . Including the shal¬ 
lower profile of MACS0329.7-0211 results in an ensemble average 
gradient of d(Z/Z 0 )(d log(r))-‘~ -0.4 (see Table[6]& Figure[l2). 

While major mergers may be the dominant mechanism by 
which the BCG grows fLidman et al.[|2013] >, we conclude that they 
are less important in the build-up of ICL outside of the BCG be¬ 
cause of the significant observed colour and metallicity gradients. 
Between 10 and 110 kpc we find total ICL luminosities of 4-8 L * 
(Table [6). If the ICL is built-up solely from stars ejected during 
accretion of L* galaxies by the BCG, and -20% of the stars of 
the merging galaxy are deposited in the ICL (as in |Murante et al.| 
( |2007| ); |Contini et ak] { |2013| )) then 20-40 merging events must oc¬ 
cur. This number of merging events is 10 times greater than that 
found in models of BCG evolution since z =1 jLidman et al.[|2013| . 
Further, this number of major mergers is enough to significantly 
erode any underlying stellar population gradient |Kobayash?]|2004[ 
Di Matteo et al. 2009). Additionally, the higher central concentra¬ 
tion of ICL relative to the galaxy distribution is consistent with a 
formation mechanism that is more efficient deeper in the cluster 
potential well e.g. tidal interactions < |Zibetti et al.[|2005| l. 

Further, the near-solar metallicities imply that the ICL cannot 
be formed solely from low-metallicity dwarf galaxies. This result 
echoes the conclusions of |Contini et al.| ( [?013| >, whose models in¬ 
dicate that galaxies with stellar masses of -10 11 M© are the primary 
contributors to the ICL. They frame their conclusion in terms of 
dynamical friction: massive satellite galaxies experience the great¬ 
est dynamical friction, causing them to quickly fall into the central 
region of the cluster where they can be tidally stripped and dis¬ 
rupted. If the observed 4-8 L* of ICL luminosity formed solely 
from dwarf disruption, then hundreds of dwarfs would be required. 
The faint end of the luminosity function must have evolved from a 
slope at least as steep as a = -1.8 to the flatter values observed to¬ 
day (-0.84 |Lm et al.|j2004| ) in order for the difference in integrated 


luminosity between the past and present luminosity functions to be 
greater than 4L*, the minimum observed ICL luminosity. However, 
the faint-end slope has been shown to have little or no evolution 
since z-1.3 (|Mancone et al.[|2012| > and possibly since z=3.5 ( |Ste-| 
|fanon & Marchesini||2013| >. 


Zibetti et al. 


( |2005| ) and |Krick et aL] < |2006| ) find that the op¬ 


tical colours of the ICL are well-matched by the cluster galaxy 
colour distribution, suggesting that the ICL results in part from stars 
stripped and disrupted from red sequence galaxies. The observed 
F110W-F160W colour range of our clusters also closely matches 
their red sequence colours, in general spanning colours equal to the 
cluster’s red sequence down to M*+2 mag. This further supports 
the idea that the progenitor population of the ICL is largely from 
L* galaxies. 

While it is possible, though unlikely, for a colour and metallic¬ 
ity gradient to develop via in situ formation, the metallicity of the 
resulting intracluster stars would likely be similar to the ICM metal¬ 
licity of ~0.3 Z 0 ( [Edge & Stewart[p^T||Sivanandam et aL]|2009| ). 
In two of our clusters (MACS 1149.6+2223 and MACS2129.4- 
0741), the metallicity gradients reach down to -0.3 Z© by 110 kpc. 
However, those metallicity profiles are produced under the assump¬ 
tion of an old stellar population with z f =3. Assuming instead that 
the ICL is developed via in situ star formation and that the pre¬ 
dicted ICL formation epoch of z<l is correct ( [Murante et al.[|2007| 
De Lucia & Blaizot 2007), the formation redshift of the intracluster 
stars would be around z f =l. For example, a stellar population at the 
redshift of MACS 1149.6+2223 (z=0.544) that formed at z=l with 
a metallicity of Z=0.4 Z© would have a FI 10W-F160W colour of 
0.21, as found using a BC03 SSP model with a burst of star forma¬ 
tion with t= 1 Gyr e-folding time and a Chabrier IMF. This is incon¬ 
sistent with the observed F110W-F160W colour profiles, where, 
for example, MACS 1149.6+2223 has a F110W-F160W colour of 
0.33 at r=110 kpc. If the ICL was formed in large part by in situ 
star formation, we would expect to see much bluer colour profiles. 


8 CONCLUSIONS 

In this paper we have presented initial results from a comprehensive 
HST program to quantify the properties of intracluster light in a 
statistical sample of clusters and groups at intermediate redshift. 
Using the spectacular spatial resolution and depth of HST, we have 
observed surface brightness and colour profiles for a unique sample 
of four galaxy clusters at redshifts 0.44<z<0.57, ranging in mass 
from 0.6 to 2.6 xl0 15 M© . 

(1) The F160W surface brightness profiles of our sample have 
an average RMS scatter between 10<r<110 kpc of 0.24 mag 
arcsec _2 (See Figure [8j. This scatter is physical - it exceeds 
the observational errors, straightforward expectations from the 
range of cluster masses in our sample, and predictions based on 
published evolutionary models for the variance attributable to 
the redshift span of our sample. We associate the additional 
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scatter with differences in ICL assembly process, formation 
epoch, and/or ICL content. 

(2) Three out of four of our clusters have intracluster stellar pop¬ 
ulations that get progressively bluer with increasing radius, 
with a mean colour gradient of -0.085+0.021 mag log(kpc) -1 
(-0.064+0.021 mag log(kpc) -1 including the flat profile of 
MACS0329.7-0211). If we interpret this colour gradient as a 
metallicity gradient for an old stellar population, the metal- 
licity decreases from high values (> Z Q ) in the central region 
where the brightest cluster galaxy (BCG) dominates (< 10 kpc) 
to subsolar at the largest radii. The average metallicity gradient 
for the 3 clusters with negative colour gradients is d(Z/Z Q )(d 
log(r)) _1 =-0.57+0.13 log(kpc)" 1 (-0.43+0.13 log(kpc)- 1 for all 
4 clusters). Such negative colour gradients can arise from tidal 
stripping of L* galaxies and/or the disruption of dwarf galax¬ 
ies, but not major mergers with the BCG. 

(3) The ICL at ~110 kpc has a colour comparable to M*+2 red 
sequence galaxies, suggesting that out to 110 kpc the ICL is 
dominated by stars liberated from galaxies with L > 0.16L*. 
The colours (and equivalently high inferred metallicities) dis¬ 
favour disruption of even lower luminosity dwarf galaxies as 
the dominant formation mechanism for the ICL. 

(4) For 10<r<110 kpc, the luminosity of the ICL for these clus¬ 
ters is 4-8 L*. For dwarf disruption to be the sole source of 
ICL, more substantial evolution of the faint end of the lumi¬ 
nosity function than observed would be required. Thus, the to¬ 
tal luminosity of the ICL also argues against dwarf disruption 
as the dominant progenitor population for the ICL. Further, the 
observed ICL luminosity also disfavours major mergers with 
the BCG as the dominant formation mechanism of the ICL as a 
merger model origin for the ICL requires 10 times more merg¬ 
ers than predicted in BCG formation models (|Lidman et al.| 

j20T2). 

To summarize, using this sample of CLASH clusters we set 
empirical limits on the amount of scatter in ICL surface brightness 
profiles at z=0.5 for theoretical models and constrain the progen¬ 
itor population and formation mechanism of the ICL. The results 
from the absolute colour, colour gradients, and total luminosity of 
the ICL, as enumerated above, are suggestive of a formation his¬ 
tory for these clusters in which the ICL is built-up by the stripping 
of >0.2L* galaxies, and disfavour significant contribution to the 
ICL by dwarf disruption or major mergers with the BCG. Our fu¬ 
ture work will include analysis of the remaining clusters from the 
CLASH sample, as well as eight groups from our HST-GO pro¬ 
gram (Program #12575), expanding the mass range of the cluster 
environments. This larger mass range, combined with the more- 
extended redshift baseline, will allow us to trace the build-up of 
the ICL over the past 5 Gyr. In addition, we will use this high- fi¬ 
delity, multi-band photometry to structurally disentangle the BCG 
and ICL. 
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